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ADMISSIBLE  TRANSLATES  OF  STABLE  PROCESSES: 


1.  INTRODUCTION 


Throughout  X  =  {X(t,w),  -oo  <  t  <  00}  is  a  stochastic  process  defined  on  a  probability  space 
and  s  =  {s(t),  -00  <  t  <  00}  is  a  nonrandom  function.  The  map  u>  — *  X(-  ,  u>)  from  the 
probability  space  to  the  space  of  all  functions  on  the  read  line,  induces  a  probability  measure  on  the 
(r-field  of  cylinder  sets,  the  distribution  of  the  process  X.  If  the  distribution  of  the  translate  of  X  by  s, 
Pg_^,  is  absolutely  continuous  with  respect  to  the  distribution  of  X,  P^,  then  s  is  called  an  admissible 
translate  of  X,  and  if  Pg_^  an<^  ^X  ate  s^n8u^ati  then  s  is  called  a  singular  translate.  The  set  of  all 
admissible  translates  of  X  is  denoted  by  AT(X). 

In  signal  detection  s  is  the  signal,  X  is  the  random  noise,  and  based  on  a  sample  observation  (a 
function  on  the  real  line)  one  decides  whether  the  observation  is  due  to  signal  plus  noise  (s+X)  or  to 
noise  alone  (X).  For  a  singular  translate  s,  in  principle,  a  correct  decision  can  be  made  with  pro¬ 
bability  one,  i.e.  signals  that  are  singular  translates  of  the  noise  can  be  perfectly  detected.  For  an 
admissible  translate  s  the  decision  is  based  on  the  Radon-Nikodym  derivative  of  Pg+X  w't^  respect  to 
Px  (likelihood  ratio)  according  to  the  Neyman-Pearson  rule  and  the  resulting  probability  of  detecting 
correctly  the  presence  of  the  signal  is  always  strictly  less  than  one.  Most  real  life  signal  detection 
problems  correspond  to  admissible  translates. 

In  Section  2  we  describe  the  well-known  complete  results  on  the  admissible  and  singular  translates 
of  Gaussian  processes.  In  Section  3  we  summarize  some  recent  results  for  stable  processes  which  are  far 
from  being  complete.  We  refer  to  [1]  and  references  therein. 

It  turns  out  that  the  widely  used  moving  average  and  Fourier  transform  models  typically  have  all 


□  □ 


2.  GAUSSIAN  PROCESSES 


Let  X  be  a  Gaussian  process  with  mean  zero  and  covariance  function  R.  The  set  of  admissible 
translates  of  X  is  precisely  the  reproducing  kernel  Hilbert  space  of  its  kernel  R,  and  all  other  translates 
are  singular.  More  concrete  representations  follow. 

Every  Gaussian  process  X  can  be  represented  in  terms  of  a  Gaussian  pro. -ss  £  with  independent 
increments  as  follows: 


(1) 


X(t)  =  f°°  f(t,A)  d£(A),  -oo  <  t  <  oo, 
J-oo 


where  f(t,  •  )  f  L2(ji)  and  ft  is  the  control  measure  of  £  :  dft(X)  =  E  |d((A)|2.  The  integral  J  g  d£ 
is  defined  for  all  g  e  L2(/j)  and  is  a  normal  random  variable  with  characteristic  function 

E  exp{  ir  |  g  d{  }  =  exp{  -  |  r2  |  |  g  |2  dft  }. 

roo  _ 

The  covariance  function  of  X  is  represented  as  R(t,s)  =  f(t, A)  f  (s,A)  dji(A).  Then 

J-oo 

foo 

(2)  AT(X)  =  {  s(t)  =  f(t,A)  g(A)  d^A),  g  e  L2(p)  }. 

J-oo 


Stationary  Gaussian  processes  (that  are  continuous  in  probability)  have  a  spectral  representation: 

(3)  X(t)  =  f°°  eiU  d£(A), 

J-oo 

where  £  has  independent  increments  and  finite  spectral  (control)  measure  p,  and  their  covariance 
roo  i(t— s)A 

function  is  R(t,s)  =  e'  d/j(A).  Their  admissible  translates  are  Fourier  transforms  of 

J-oo 

finite  signed  measures  that  are  absolutely  continuous  with  respect  to  ft  with  p-square  integrate  Radon- 
Nikodym  derivative: 

(4)  AT(X)  =  {  s(t)  =  eitA  di/(A),  u«fi,  ^  c  L 2(ft)  }. 

When  ft  <<  Leb  with  spectral  density  ^(A)  =  d^(A)/dA  then 

AT(X)  =  {  s(t)  =  eitA  S( A)  dA,  dA  <  oo  }, 
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and  in  addition  to  the  spectral  representation  (3),  X  also  has  a  moving  average  representation: 

roo 

(5)  X(t)  =  h(t-s)  d*(s), 

J-oo 

where  in  (5),  £  has  stationary  independent  increments  (i.e.,  n  =  Leb)  and  h  e  L2;  its  covariance 

function  is  R(t,s)  =  [  e*^  |h( A)|2  dA  ;  and  its  admissible  translates  are 

J  -co 

roo 

(6)  AT(X)  =  {  s(t)  =  h(t-s)  g(s)  ds,  g  f  L2  }. 

J  -oo 

For  nonstationary  processes  representations  in  terms  of  processes  with  dependent  increments  are 
useful.  Every  Gaussian  process  X  can  be  represented  in  terms  of  a  Gaussian  process  tj  with  possibly 
dependent  increments  as  follows: 

roo 

(l')  X(t)  =  f(t,u)  drj(u),  -oo  <  t  <  oo. 

J-oo 

Here  tj  is  an  L2(fi,<J,P)-valued  Gaussian  measure  generally  not  independently  scattered,  with  bimeasure 
0:  0(B,C )  =  E{  r?(B)  rj( C)  },  B,C  :  bounded  Borel  sets,  and  each  f(t,  •  )  is  77-integrable.  The  co- 
variance  of  X  is  R(t,s)  =  ||  f(t,u)  f  (s,v)  d/3(u,v),  and  its  admissible  translates  are 

.  f  f°°  f°° 

(2')  AT(X)  =  {  s(t)  =  j  J  ^  f(l-u)  «(v)  d0(u>v)  =  j  ^  f(l<u)  d/*g(u)-  8  :  »?-integrable  } 

roo 

where  /3g  is  the  measure  0g(-)  =  J  g(v)  0( •  ,dv). 

When  f(t,u)  =  eltu  in  (]/)  and  rj  has  bounded  semi- variation,  then  X  is  harmonizable: 


with  bispectral  measure  0  and  covariance 


roo 

X(t)  =  eltu  dr?(u), 

J-oo 

_.  ,  f  f00  i(tu— sv) 

nee  R(t,s)  =  JJ  e 


d/?(u,v).  It  is  stationary  only  if 


its  bispectral  measure  0  is  concentrated  on  the  diagonal  of  the  plane.  Its  admissible  translates  are: 

(4f)  AT(X)  =  {  s(t)  =  [  eltu  d/?R(u),  g  :  q-integrable  }. 

J-00  0 

When  f(t,u)  =  h(t— u)  in  (l*),  X  is  a  “generalized”  moving  average  (of  a  dependent  increments 
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process): 

,  roc 

(5  )  X(t)  =  h(t-u)  dr;(u), 

J-oo 

whose  admissible  translates  are 

l  foo 

(6;)  AT(X)  =  {  s(t)  =  J  h(t-u)  d/?g(u),  g  :  q-integrable  }. 

The  representations  of  the  admissible  translates  of  general  nonstationary  processes,  (2*),  (4*),  and 
(61),  are  similar  to  those  of  stationary  processe,  (2),  (4)  and  (6),  but  not  as  explicit  and  simple  because 
generally  no  explicit  description  of  the  rj-integrable  functions  is  available. 

3.  STABLE  PROCESSES 

In  this  section  X  is  a  symmetric  a-stable  (SaS)  process  (0  <  a  <  2),  i.e.  all  linear  combinations 
c1X(t1)+  •••  +cnX(tn)  are  SaS  random  variables  (with  characteristic  functions  of  the  form 
exp{  -  const.  |  r  |a  }  ).  When  a  =  2,  X  is  Gaussian.  When  0  <  a  <  2,  X  is  non-Gaussian  stable 
and  can  be  represented  as  in  (1),  where  £  has  SaS  independent  increments  and  control  measure  fi. 
Here  the  integral  J  g  d£  is  defined  for  all  g  e  La(/i)  and  is  a  SaS  random  variable  with 
characteristic  function 

E  exp{  ir  J  g  d£  }  =  exp{  -|  r  |a  j  |  g  |a  dp  }. 

A  crucial  difference  in  the  representation  (1)  between  the  Gaussian  (a=2)  and  non-Gaussian 
(0  <  a  <  2)  cases  is  that  whereas  in  the  Gaussian  case  the  linear  space  of  the  increments  of  £  can 
always  be  taken  to  be  equal  to  the  linear  space  of  the  process  X,  in  the  non-Gaussian  case  the  former 
is  generally  larger  than  the  latter  (and  generally  infinite  dimensional  even  when  the  linear  space  of  X  is 
finite  dimensional). 

The  analog  of  (2)  is  no  longer  valid  when  0  <  a  <  2  and  all  that  can  be  said  in  general  is  that 

(7)  AT(X)  c  {  s(t)  =  f°°  f(t,A)  g(A)  d^(A),  g  e  L  ,(»)  }  £  F(X), 

J-oo 
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where  1/a  +  1/a*  =  1,  and  that  translates  outside  the  set  F(X)  are  singular.  AT(X)  may  be  as 
large  as  F(X)  (e.g.  when  X  is  sub-Gaussian)  or  as  small  as  {0}  (e.g.  when  X  is  Levy  motion). 

In  order  to  describe  better  the  contrasts  and  similarities  between  the  Gaussian  and  other  stable 
cases  we  first  concentrate  on  the  independently  scattered  SaS  measure  £  =  {  £(B),  B  e  31^}  in  (1), 
where  denotes  the  bounded  Borel  sets  of  the  real  line.  When  a  =  2,  AT(£)  consists  of  all  signed 
measures  which  are  absolutely  continuous  with  respect  to  p  with  Radon-Nikodym  derivative  in  L2(p). 
In  sharp  contrast  when  0  <  a  *  2  and  p  is  nonatomic,  then  £  has  no  admissible  translates  and  in  fact 
all  translates  are  singular.  On  the  other  hand,  when  0  <  a  <  2  and  p  is  purely  atomic  with  atoms 
{  an}„—  p  1  <  N  <  oo,  then  AT(£)  consists  of  all  measures  s  concentrated  on  { an } ^  with 

fi  N{an})|2  /  p2/“({an})  <oo; 
n=l 

this  result  is  due  to  Shepp  [  2  ]  and  is  valid  for  all  0  <  a  <  2.  For  a  general  control  measure  AT(£) 

=  {0}  when  p  has  no  atoms,  and  when  fi  has  atoms  {an}^_  j  then  AT(£)  is  as  above. 

We  now  consider  a  general  SaS  process  with  representation  (1)  and  0  <  a  <  2.  X  is  called 
invertible  if  its  representors  {f(t,  •  ),  -oo  <  t  <  oo  }  are  complete  in  La(p),  i.e.  if  the  linear  spaces 
of  X  and  of  the  increments  of  £  are  equal.  When  X  is  invertible  then  every  translate  is  either  ad¬ 
missible  or  singular  and  AT(X)  =  {0}  when  p  has  no  atoms  and  if  \i  has  atoms  {an}^_|,  1  <  N  <  oo, 
then 

N  1/a  N 

AT(X)  =  {  s(t)  =  E  s n  n  '  ({an})  f(t,an)  ,  £  |  sn|2  <  oo  }  . 

n=l  n=l 

Here  AT(X)  is  a  projier  subset  of  F(X)  even  when  fi  is  purely  atomic. 

When  X  is  stationary  and  has  the  spectral  representation  (3),  where  £  has  independent  SaS 

increments  and  finite  spectral  (control)  measure  p,  then  X  is  invertible.  In  the  Gaussian  case  its 

admissible  translates  are  described  in  (4).  In  the  non-Gaussian  case  0  <  a  <  2,  if  the  spectral  measure 

N 

p  has  no  atoms  all  translates  of  X  are  singular  and  if  p  has  atoms  {an}j,_.i>  1  <  N  <  oo,  then 

AT(X)  =  {  s(t)  =  fi  sn  p1/Q({an})  eiU"  ,  fi  |  sn|2  <  oo  }. 
n=l  n=l 
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When  X  is  stationary  and  has  a  moving  average  representation  (5),  where  £  has  stationary 
independent  SaS  increments,  then  all  translates  of  X  are  singular  when  0  <  a  <  2,  provided  the 
translates  of  its  kernel  h  are  complete  in  La.  In  the  Gaussian  case  the  admissible  translates  of  X  are 
described  in  (6). 

The  processes  with  spectral  (3)  and  those  with  moving  average  (5)  representation  are  the  most 
widely  studied  classes  of  stationary  stable  processes,  and  are  actually  disjoint  when  0  <  a  <  2,  while 
in  the  Gaussian  case  a  =  2  the  latter  is  a  subset  of  the  former.  Since  when  0  <  a  <  2  in  most  cases 
they  have  no  admissible  translates  and  all  their  translates  are  singular,  they  do  not  provide  reasonable 
models  of  stable  noise  in  signal  detection,  as  every  signal  can  be  detected  perfectly  in  principle  in  the 
presence  of  such  noise.  One  way  to  introduce  more  realistic  stable  noise  models  is  to  allow  the 
increments  of  £  in  (3)  and  (5)  to  have  some  dependence,  i.e.  to  consider  the  stable  counterparts  of  (3,J 
and  (§'). 

Thus  let  X  have  representation  (1,)1  where  ij  has  dependent  SaS  increments  with  0  <  a  <  2. 
Here  rj  is  a  SaS  measure  which  is  not  independently  scattered  and  can  be  represented  as 

7(-)=  f°°  "('  ,A)d{(A), 

J-oo 

where  (  has  independent  SaS  increments  and  control  measure  /j,  and  v  is  an  La(/r)-valued  measure.  If 

g  is  j/-integrable  then  J  g  d v  t  La(/i)  and 

roo  roo  roo 

g  d»?  =  {  g  di/  }  d{. 

J-OO  J-oo  J-oo 

Thus  with  each  f(t,  •  )  being  i/-integrable  we  have 
roo  roo  roo 

(8)  X(t)  =  f(t,u)  d??(u)  =  {  f(t,u)  t/(du,A)  }  d£(A),  -oo  <  t  <  oo. 

J-oo  J-oo  J-oo 

To  ensure  that  r?  is  not  independently  scattered  we  assume  that  for  some  disjoint  Borel  sets 

and  B2,  fi{  A  :  j^B^A)  i/( B2,A)  ^  0  }  >  0,  since  »7(Bj)  and  »?(B2)  are  independent  iff  t/(Bj,A)  i/(B2,A) 

roo 

=  0  a.e.  hi].  Note  that  »?({a})  =  0  iff  |//({a},u)|2  d/i(u)  =  0,  so  that  t]  may  have  no  atoms  even 

J-oo 

when  £  does,  e.g.  if  i>  has  no  atoms:  i/({a},u)  =  0  a.e.  [p].  It  is  therefore  possible  for  X  to  have 
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admissible  translates  even  when  rj  has  no  atoms. 


roo 

We  further  require  X  to  be  invertible,  i.e.  {  f(t,u)  ^(du,  •  ),  -oo  <  t  <  00}  to  be  complete 

J-00 

N 

in  La(/i),  and  /r  to  have  atoms  {an}^_p  1  <  N  <  00.  It  then  follows  that 

N  i/n  foo  N  , 

AT(X)  =  {  s(t)  =  £  sn  n  '  ({an})  f(t,u)  i/(du,an)  ,  £  lsnl  <  0°  }• 

n=l  J'00  n=l 

.  l’°° 

When  N  <  00  the  representation  of  an  admissible  translate  can  be  written  as  in  (2  ):  s(t)  =  f(t,u) 

J  -00 


di/jy(u)  where  the  measure  is 

(9) 


•  )  =  L  sn  ^'^({an})  H  '  .  an)- 
n=l 


When  N  =  00  the  infinite  series  in  (9)  does  not  generally  define  a  measure  and  we  can  only  write  s(t)  = 

f°° 

lim^^^l  f(t,u)  di/j^u),  where  1/^  is  defined  as  in  (9). 

With  f(t,u)  =  eltu  in  (8)  and  r)  of  bounded  semi- variation  then  X  is  harmonizable: 

X(t)  =  [  eltu  dij(u)  =  [  d£(A), 

J-00  J-00 

where  £>(t, A)  =  [  e*tu  i/(du,A),  with  admissible  translates  all  functions  of  the  form  (as  in  (4')) 

J-00 

roo  ..  roo  . 

s(t)  =  eltu  di/N(u)  if  N  <  oo,  =  lim  eltu  di/.  (u)  if  N  = 

J-00  1  k— 1 -oo  J-00 


00  . 


Harmonizable  processes  are  nonstationary  when  r)  has  dependent  increments. 

When  f(t,u)  =  h(t  — u)  in  (8),  X  is  a  “generalized”  moving  average  (of  a  dependent  increments 
process): 

roo  roo  roo 

(10)  X(t)  =  h(t-u)  drj(u)  =  {  h(t-u)  t/(du,A)  }  d£(A) 

J-00  J-00  J-00 

with  admissible  translates  are  ail  functions  of  the  form  (as  in  (6*)) 

roo  foo 

(11)  s(t)  =  h(t-u)  dirN(u)  if  N  <  00,  =  lim  h(t-u)  di/,(u)  if  N  =  00. 

J-00  k  — ►  00  J-00 

It  is  easy  to  construct  examples  of  nonstationary  generalized  moving  averages,  but  there  exist  sta- 
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tionary  ones  as  well.  Indeed,  taking 


i/(B„\)  =  0(A)  f  e‘x9(A)  dx 
JB 

roo 

for  bounded  Borel  sets  B  where  4>  c  La(/i)  we  have  for  g  e  L1?  g(x)  i/(dx,A)  =  <^(A)  g[0(A)] 

J  -oo 

and  from  (10), 


X(t)  =  [°°  e'te(A)  <t>(\)  h[-tf(A)]  d((A). 
J-oo 


Therefore  X  is  stationary  and  invertible  with  admissible  translates  all  functions  of  the  form  (as  in  (6) 
and  (6*)) 


(t)  =  (  h(t-u)  g(u)  du, 
J  -oo 


where  g(u)  =  £  sn  ^“({a,,})  <6(an)  eiu^an)  and  £  |sn|2  <  oo  . 
n=l  1 

Thus  harmonizable  and  generalized  moving  averaged  SaS  processes  can  be  constructed  with 
rich  classes  of  admissible  translates  which  may  serve  as  realistic  models  of  noise  in  signal  detection. 
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